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Abstract 
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processes involving rj and rj' include a large flavor-singlet penguin amplitude (5) as proposed earlier 
and a penguin amplitude Ptu associated with intermediate t and u quarks. For the B^ — > tt+t/' 
mode a term Stu associated with intermediate t and u quarks also may be needed. Values of the 
weak phase 7 are obtained consistent with an earlier analysis of ^ VP decays, where V denotes 
a vector meson, and with other analyses of CKM parameters. 
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I. INTRODUCTION 



A central objective of the study of B meson decays is to help determine the phases and 
magnitudes of Cabibbo-Kobayashi-Maskawa (CKM) matrix elements, through the measure- 
ment of branching ratios and CP-violating observables. It is important to have accurate 
and self-consistent information on CKM matrix elements if they are ever to be compared 
with fundamental theories predicting them. At present no such theories exist. A further 
objective is to learn about possible new physics at higher mass scales, affecting rare B de- 
cays by giving observables that appear to be inconsistent with others. One wishes to know 
whether there are any sources of CP violation other than the phases in the CKM matrix 
first proposed by Kobayashi and Maskawa 

Charmless B meson decays, many of whose branching ratios and CP asymmetries (CPA's) 
have been measured to good accuracy, are an interesting and useful set of modes. Following 
the method presented in Ref. P for B decays into a vector meson {V) and a pseudoscalar 
meson (P), we analyze observables in B decays into two pseudoscalar mesons {B PP 
decays) in the present paper. From the results of fits involving a small set of invariant 
amplitudes, one can extract information about the parameters in theory, compare with other 
known constraints, and predict as-yet-unreported observables. In particular, the amplitudes 
contributing to two-body hadronic charmless B decays involve only one nontrivial weak 
phase 7 within the standard model (SM). In a previous analysis of P — > VP decays P|, 
we found good agreement between the favored range of 7 from a fit to the VP modes and 
that from fits to CKM parameters ^| based on other measurements. It is therefore of great 
interest to see if the PP modes give a consistent result. 

In p.e.ent analy... we flavo. SU,3) .y,„n,e.. QBBBBQQ a. a wo.W 
hypothesis. Motivated by factorization in tree-level amplitudes, we take symmetry breaking 
due to decay constant differences into account in these amplitudes when relating strangeness- 
conserving and strange-changing processes. We leave the issue of SU(3) symmetry break- 
ing in penguin- type amplitudes to experimental data. As a test, one can compare the 
P+ 71^ mode (involving purely a strangeness-changing QCD penguin amplitude) with 
the B^ K^K^ and P° — K^K^ modes (involving purely strangeness-conserving QCD 
penguin amplitudes). In the limit of flavor SU{'i) symmetry, they should differ by a ratio 
of CKM factors, Vcs/Vcd- If penguin amplitudes Ptu associated with intermediate t and u 
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quarks are important, the predictions for these modes will be affected. 

We find acceptable fits to B ^ mi and B nK branching ratios and CP asymmetries 
with a combination of tree, color-suppressed (C), penguin (P), and electroweak penguin 
{Pew) amplitudes. In contrast to an earlier analysis of B ^ PP decays in order to 
describe these decays we must introduce a rather large value of \C/T\ and a non-trivial 
relative phase between C and T. A large \C/T\ value could improve agreement between the 



12| . Our conclusion is driven in part by the 
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QCD factorization approach and experiment 
large branching ratio for B^ 7r°7r° reported recently 

The data on processes involving rj and t]' also have made some progress since our earlier 
analysis Crucial additional terms for describing these decays include not only a large 
flavor- singlet penguin amplitude (S) as proposed (e.g.) in Refs. but also a penguin 
amplitude Pm associated with intermediate t and u quarks, and (for the B^ —>■ n^Vj mode) 
a term St^ associated with intermediate t and u quarks. 

Values of the weak phase 7 ~ 60° are obtained consistent with our earlier analysis of 
B VP decays [2] . Other robust aspects of our fit include the magnitude of the strangeness- 
changing penguin amplitude, the strong phase of the tree amplitude relative to the penguin 
(~ 20°-30°), the size of electroweak penguin contributions, the correlation of a large direct 
CP asymmetry in B^ tt+tt^ with a small one in B^ 7[~K~^, the correct prediction of 
signs and magnitudes of all other measured direct CP asymmetries as well, and a fairly large 
negative value of the time-dependent CP asymmetry parameter St,-,,. Some other aspects of 
the fit are less likely to remain unchanged in the face of further data; we shall comment on 
them in due course. 

We review our conventions for the quark content of pseudoscalar mesons and topological 
amplitudes in Section |nj Experimental data and topological decompositions of decay ampli- 
tudes are presented in Section UTTl In Section |TVl we enumerate the data that will be used in 
our fit- Two fits to tttt and ttK observables are presented in Section |3 while modes with 
1] or 1]' in the final state are included in Section IVII We comment on robust and less-stable 
aspects of the fits in Section IVTTl Based upon our fitting results, we discuss our predictions 
for as-yet-unreported modes in Section IVIIII Comparisons with other recent approaches 
(e.g., Refs. QMHQ) are pursued in Section HXl We summarize our findings in Section 

m 
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II. NOTATION 



Our quark content and phase conventions T], Isf are: 

• Bottom mesons: = db, B^ = bd, B^ = ub, B^ = —bu, Bg = sb, Bg = bs; 



Charmed mesons: = —cu, = uc, = cd, D = dc, Dt^ = cs, D = 



sc: 



• Pseudoscalar mesons: = ud, = {dd — uu)/\^, vr = —du, = us, = ds, 
= sd, = —su, T] = {ss — uu — dd) / -\/3, rj' = {uu + dd + 2ss) / -\/6; 

The rj and 1]' correspond to octet-singlet mixtures 

Tj = r]s cos - r]i sin Oq , r]' = r]s sin + r]i cos 6*0 , (1) 

with ^0 = sin-i(l/3) = 19.5°. 

In the present approximation there are seven types of amphtudes: a "tree" contribution 
t; a "color-suppressed" contribution c; a "penguin" contribution p; a "singlet penguin" 
contribution s, in which a color-singlet qq pair produced by two or more gluons or by a Z or 
7 forms an SU{3) singlet state; an "exchange" contribution e, an "annihilation" contribution 
a, and a "penguin annihilation" contribution pa. These amplitudes contain both the leading- 
order and electroweak penguin contributions, and appear in the independent combinations 

t = T + Pgw ! c = C + Pew , 

p = P — Ptu — |Pew ; s = S — Stu — |Pew , (2) 
a = A , e + pa = E + PA , 

where the capital letters denote the leading-order contributions , Q while Pew and 



P^^JSl are respectively color-favored anc 
8|. We shall neglect smaller terms 
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color-suppressed electroweak penguin amplitudes 
PE-Ew and PA-^w [{l, Z)-exchange and (7, Z)- 
direct-channel electroweak penguin amplitudes]. We shall denote AS* = transitions by 
unprimed quantities and \^S\ = 1 transitions by primed quantities. The hierarchy of these 
amplitudes can be found in Ref. 21[. By writing QCD and flavor-singlet penguins as P — Pm 
and S — Stu, we adopt the so-called c-quark convention, in which the heavy top quark is 
integrated out from the theory. For penguin-type amplitudes, we use the unitarity relation 
Vtb^td + Vcb^cd + Vub^ud = to remove any top quark dependence. The V^*feV^d term of the 
top quark mediated penguin is combined with the up quark mediated penguin to form Ptu or 



Stu- Similarly, the V*^cd term is united with the charm quark mediated penguin into P or 
S. As a consequence, the strangeness-conserving P and 5* and strangeness-changing P' and 
S' penguin amplitudes have real weak phases in our discussions. The relation between the 
c-quark convention and the t-quark_convention, where the c quark dependence is removed 
instead, can be found in, e.g., Ref. 

The partial decay width of two-body B decays is 



V{B M1M2 



Pc 



■\A{B M1M2] 



(3) 



STrm^ 

where Pc is the momentum of the final state meson in the rest frame of B, itlb is the B 
meson mass, and Mi and M2 can be either pseudoscalar or vector mesons. Using Eq. Q, one 
can extract the invariant amplitude of each decay mode from its experimentally measured 
branching ratio. To relate partial widths to branching ratios, we use the world-average 
lifetimes r+ = (1.653 ± 0.014) ps and = (1.534 ± 0.013) ps computed by the LEPBOSC 
group [2^. Unless otherwise indicated, for each branching ratio quoted we imply the average 
of a process and its CP-conjugate. 



III. AMPLITUDE DECOMPOSITIONS AND EXPERIMENTAL RATES 



The experimental branching ratios and CP asymmetries on which our 



are listed in Tab 
29 
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es m and |IT1 Contributions from the CLEO 
37|, and Belle 14, 36, 38, 3, 40, 41, 41 
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analysis is based 
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, BaBar 







46| Collaborations 



are included j42|. In order to implement upper bounds in a consistent manner we have 
computed our own experimental averages for B^ vr+r/' and B^ r]K^. These two modes 
were observed by BaBar with a significance of 3.4 and 3.3 standard deviations, respectively. 

We list theoretical predictions and averaged experimental amplitudes for charmless B 
PP decays involving AS" = transitions in Table UTTl and those involving \ AS\ = 1 transitions 
in Table Hvl Theoretical predictions are shown in terms of topological amplitudes t, c, p and s 
while e, a and pa contributions are neglected. They are expected to be suppressed by a factor 
of order 1 / rrih relative to tree and penguin amplitudes A suppression factor proportional 
to fs/^h was suggested in 0, Q]. Future measurements of the B^ K~ decay mode 

which only receives contributions from exchange and penguin annihilation diagrams will test 
this suppression. 
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TABLE I: Experimental branching ratios of selected AS = decays of B mesons. CP-averaged 
branching ratios are quoted in units of 10^^. Numbers in parentheses are upper bounds at 90 % 
c.l. References are given in square brackets. Additional lines, if any, give the CP asymmetry Acp 
(second line) or {S,A) (second and third lines) for charged or neutral modes, respectively. The 
error in the average includes the scale factor S when this number is shown in parentheses. 

Mode CLEO BaBar Belle Average 

B+ ^ TT+vr" 4.6+1-*+°;^ 124] S.slg ° ± 0.6 [2^] 5.0 ± 1.2 ± 0.5 |14] 5.2 ±0.8 

-0.03to;l? ± 0.02 f2§] -0.14± 0.24l°;°4 [38] -0.07 ±0.14 

K+K° < 3.3 [24] 1.1 ± 0.75+q {1 (< 2.5) [29] < 3.3 jl4] < 2.5 

TT+r] 1.2l^;^ (< 5.7) [25] 5.3 ± 1.0 ± 0.3 [30| 5.4t^;° ± 0.6 [39] 4.9 ± 0.9 

-0.44 ±0.18 ±0.01 [30] - -0.44 ±0.18 

TT+V l.oli;o (< 12) [25] 2.7 ± 1.2 ± 0.3 {< 4.5) [22] < 7 [4^] 2.4 ± 1.1 (< 4.5) 

^ TT+TT- 4.5+j;2lo;4 [24] 4.7 ±0.6 ±0.2^ 4.4 ± 0.6 ± 0.3 fl4] 

{-0.40 ±0.22 ±0.03 I -1.00 ±0.21 ±0.07 

0.19 ± 0.19 ± 0.05 y 0.58 ± 0.15 ± 0.07 

-n-O-n-O < 4.4 [24] 2.1 ± 0.6 ± 0.3 [13] 1.7 ±0.6 ±0.2 [14] 

K+K- < 0.8 [24] < 0.6 [31] < 0.7 [14] 

K°K° < 3.3 [24] OStol ± 0.1 {< 1.8) [29] < 1.5 [14] 

Tr°V 0.0^0,1 (< 2.9) [25] 0.7^11 ± 0.3 {< 2.5) [33] 
7r%' 0.0+J;^ (< 5.7) [25] 1.0+^;;^ ± 0.8 {< 3.7) [33] 

W < 18 ^ -0.9l^;!^ ± 0.7 (< 2.8) [24] - 

rm' < 27 [27] 0.6j;^;J, ± 1.1 {< 4.6) [34] 

jj'jj' < 47 [27] 1.713;* ±0.6 (< 10) [34] 

IV. FIT AND DATA POINTS 

We define for n experimental observables Xi ± AXj and the corresponding theoretical 
predictions X^^, 

The data points are the branching ratios and the CP asymmetries. We write the correspond- 
ing theoretical predictions in terms of topological amplitudes and extract their magnitudes, 
weak phases and strong phases by minimizing x^- 

Tables |l] and O contain a total of 26 data points, including 9 observables from AS = 
decays and 17 from \/S.S\ = 1 decays. The modes involving vrvr and ttK consist of the 
following 15 pieces of data: 

• The TT+TT^ decay involving the t and c amplitudes provides two data points. Since both 




TABLE II: Same as Tabled for \/^S\ = 1 decays of B mesons. 



Mode CLEO BaBar Belle Average 

S+ TT+K" iS-Stfltl'l IMl 22.3 ± 1.7 ± 1.1 [29] 22.0 ± 1.9 ± 1.1 [14] 21.8 ±1.4 

0.18 ± 0.24 ± 0.02 [2gJ -0.05 ± 0.08 ± 0.01 [29] 0-07to ogto 23 [42] 0.02 ± 0.06 

7tOK+ 12.9l^;^l^^-j [24] 12.8l};5 ± 1.0 [28] 12.0 ± 1.3^ [14] 12.5 ± 1.0 

-0.29 ± 0.23 ± 0.02 [26] -0.09 ± 0.09 ± 0.01 [28] 0.23 ± Clltg ^4 [38] 0.00 ± 0.12 (5 = 1.79) 

riK+ 2.2+2 J {< 6.9) [25] 3.4 ± 0.8 ± 0.2 [30] 5.3+^-5 ± 0.6 [39] 3. 7 ±0.7 

-0.52 ±0.24 ±0.01 [30] - -0.52 ± 0.24 

T]'K+ 80^1° ± 7 76.9 ± 3.5 ± 4.4 [2^ 78 ± 6 ± 9 [4^] 77.6 ± 4.6 

0.03 ±0.12 ±0.02 [26] 0.037 ± 0.045 ± 0.011 [^ -0.015 ± 0.070 ± 0.009 [44] 0.02 ± 0.04 

B° ^ 7r-A'+ 18.0l2 ito g [24] 17.9 ± 0.9 ± 0.7 [31] 18.5 ± 1.0 ± 0.7 [14] 18.2 ± 0.8 

-0.04 ± 0.16 ± 0.02 [2^ -0.107 ± 0.041 ± 0.013 [2S] -0.088 ± 0.035 ± 0.018 [^fi] -0.09 ± 0.03 

7T<^K° 12.8+3;°+^;^ [24] 11.4 ± 1.7 ±0.8 [22] 11.7 ± 2.3+i;3 [14] 11.7 ±1.4 



i+0.38 , 



0.48 ± 0.42 



[37] 

I -0.40to;2? ± 0.09 I -0.40 ±0.29 

r]K° O.Otolo {< 9.3) [2^ 2.9 ± 1.0 ± 0.2 (< 5.2) [30] < 12 [4g 2.5 ± 1.0 (5 = 1.08) (< 5.2) 

t^'RO 89+}g±9[25] 60.6 ± 5.6 ± 4.6 [35] 68 ± lO^g ^ 65.2 ± 6.2 (S = 1.03) 

f 0.02 ±0.34 ±0.03 f 0.43 ± 0.27 ± 0.05 f 0.27 ±0.21 

< [35] <^ [46] <^ 

-0.10 ± 0.22 ± 0.04 -0.01 ±0.16 ±0.04 -0.04 ±0.13 



amplitudes have the same weak phase except for a small contribution from EWP, no 
significant CPA is expected. 

• The TT+TT" decay involves the t and p amplitudes with different weak phases. Time- 
dependent CPA's have been observed by both BaBar and Belle groups. Thus, this 
mode provide three data points. 

• The 7r°7r° decay involving the c and p amplitudes only provides one data point because 
no CPA has been measured yet. 

• The TT^K^ decay involving only the p' amplitude provides two data points, although 
no significant CPA is expected. This mode plays a dominant role in constraining the 
magnitude of the P' amplitude. 

• The n'^K^ decay involving the p', t', and c' amplitudes provides two data points. 

• The 7r~K^ decay involving the p' and t' amplitudes provides two data points. 

• The TT^K^ decay involves the p' and c' amplitudes. Time-dependent CPA's have been 
reported by the BaBar group. Thus, this mode provides three data points. 
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TABLE III: Summary of predicted contributions to AS = decays of B mesons to two pseu- 
doscalars. Amplitude magnitudes |Aexp| extracted from experiments are quoted in units of eV. 

Mode Amplitudes pc (GeV) l^expT 

B+ ^ 7r+7r° ~72(* + ^) ^-^^^ 23.4 ±1.7 

K+K° p 2.593 < 16.4 

vr+r/ --^{t + c + 2p + s) 2.609 22.9 ± 2.0 

vr+V -^{t + c + 2p + is) 2.551 16.2 ± 3.8 (< 22.2) 

vr+vr- -{t+p) 2.636 22.8 ± 1.1 

vrOvrO ~72^^~P^ ^'^"^^ 14.7 ± 1.8 

K+K- -{e + pa) 2.593 < 8.3 

K'^K° p 2.592 < 13.2 

TT^r] -^{2p + s) 2.610 < 17.0 

TT^r]' ^{p + 2s) 2.551 < 20.9 

r/r/ ^{c + p + s) 2.582 < 18.1 

r/V -^(c + p+|s) 2.522 < 23.4 

r/V i72('^ + ^' + ^^) ^-^^0 <35.0 
" l^expl is defined by Eq. @ as an amplitude related to a CP-averaged branching ratio quoted in Table I. 

TABLE IV: Same as Table |llT| for \AS\ = 1 decays of B mesons. 

Mode Amplitudes pc (GeV) l^expl 

B+ vt+KO p' 2.614 48.2 ± 1.6 

Tr^K+ -^{p' + t' + c') 2.615 36.6 ± 1.5 

r]K+ -^{s' + t' + c') 2.588 19.9 ± 1.9 

r]'K+ -^(3p' + 4s' + t' + c') 2.528 92.5 ± 2.7 

TT-K+ -{p' + t') 2.615 45.7 ±1.0 

7r°i^° TI^^'"'^') ^-^^^ 36.6 ±2.2 

r/K° --l=(s'±c') 2.587 17.0 ± 3.5 (< 24.6) 

r/'K° ^(3p' ±4s' + c') 2.528 88.0 ± 4.2 
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Successful SU(3) fits to modes with an 77 or r( in the final state require amplitudes beyond 
those mandated by the vrvr and i^K fits. A common feature of these modes, for example, 
is that they involve a fiavor singlet amplitude s or s'. Moreover, uncertainties in 77 and rj 
wave functions and possible SU(3) breaking effects can affect such fits so we list these 
11 data points separately: 

• The TT"*"?] mode involving the combination t + c + 2p + s provides two data points. 

• The TT"'"?]' mode involving the combination t + c + 2p + 4s provides one data point. 

• The riK^ mode involving the combination s' + 1' + d provides two data points. Note 
that it does not contain p'; all three contributing amplitudes are comparable in size. 
One generally expects significant CPA as a result of the interference between tree-level 
and penguin-loop diagrams. 

• The riK^ mode involving the combination s' + d provides one data point. 

• The rjK^ mode involving the combination 3p' + 4s' + 1' + c' provides two data points. 

• The r(K^ mode provides three data points, including the CP-averaged branching ratio 
and time-dependent CPA's. 



V. FIT TO vrvr AND i:K MODES 

To avoid complication from uncertainties in the fiavor-singlet amplitudes, wave functions 
of ?7 and 77', and associated SU(3) breaking effects, we first fit the fifteen vrvr and vri^ data 
points. A study restricted to 5 ^ Kn decays based on similar assumptions was carried out 
in Refs. |4^. Guided by the relative importance of strangeness-conserving and strangeness- 
changing transitions, we choose T, C, P', and Pm as our parameters. 

We further fix the strong phase convention to be 

T = |T| 6^^'^^+^) , (5) 

C = \C\ d^^T+^C-rl) ^ (g) 

Pin = \Ptu\d^'^^^^^^ , (7) 
P' = -|P'| . (8) 
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The phase convention is such that zero strong phases of T, C and Pm amphtudes correspond 
to these amphtudes having a phase of 7 with respect to the penguin-type amphtude P. Note 
that 5c is defined as a relative strong phase between the C and T amphtudes. The extra 
minus signs for P' comes from the relative weak phase vr between P' = (Ks/Kd) P and P 
amplitudes. 

The expressions for the T', C, P/„, P, and Pew are obtained from the above equations 
taking into account the following ratios 

T' Vus fx A fx 



T Vud U 1 - A72 U 

C ^ Pla ^ Ks ^ ^ 

C ~ Ptu~ Vud~ I- AV2 
P Pew Vcd A 



0.281 , (9) 
0.230 , (10) 
-0.230 , (11) 



P' P'e^ 1-AV2" 
where A = 0.224 [50|. Therefore, a major SU(3) breaking effect from the decay constant 
difference is included for tree-type diagrams. No such effect is considered for penguin- 
type amplitudes because we do not expect factorization to work in such cases. The ratio 
Pew/P'ew — '^cd/Vcs is being used for the simplicity of our analysis. We checked that using 
Pew / P'ew — ^td/Vts (to express t-quark dominance of EWP amplitudes) does not affect the 
results in any significant way. 

We explore two approaches to fitting tttt and nK data points. One of them (Fit II) uses 
Eqs. (j21) for the topological amplitudes t, c and p: 

t^T + P^^, (12) 

c = C + Pew, (13) 

p = P-Pi„-^P^w • (14) 

Using these three equations, we can write the amplitude for any vrvr or ttK decay mode in 
Tables UTTl and HVl in terms of 9 parameters: weak phase 7, topological amplitudes |T|, \C\, 
\Ptu\, and |P'|, strong phases 6t, Sq, and 6p^^, and a parameter Sew- The latter relates EW 
penguins to tree-level diagrams and will be defined below. 

In the other approach (Fit I) we use the fact that Pm has the same weak factors as tree- 
level amplitudes T and C. This allows us to absorb the Ptu penguin into redefined T and C 
amplitudes: 

f = T-Ptu, (15) 
10 



C = C + Ptu 



(16) 



By writing topological amplitudes t, c and p in terms of T, C and P as 

t = T + PEw, (17) 
c = C + Pew , (18) 
p = P - iplw , (19) 

we still get the correct expressions for P ^ tttt and B 7iK decay amplitudes, except for 
P+ 7r+/i° and P+ tt^K~^. In these two cases the resulting expressions differ from 
the correct ones by a P/^ term. Compared to the dominant QCD penguin P', this term is 
expected to be small. Thus, Fit I gives a good description of P ^ tttt and P tcK modes 
in terms of redefined tree-level amplitudes. The advantage of this approach is a smaller 
number of fit parameters as both \Ptu\ and its strong phase Sp^^ are absorbed into T and C. 
Just 7 fit parameters are used in Fit I: weak phase 7, amplitudes |T|, \C\, and |P'|, strong 
phases 6f and 5^, and the 6ew parameter. 

A relation between the EWP amplitudes and the tree-type diagrams has been found 
in Ref. 0] using Fierz transformation to relate EWP operators with tree-level operators. 
Explicitly, we have the relations 

P'ew = -^Ew T'e-'^ = -5ew \T'\ e''^ , (20) 

P'ew = -Sew C'e-'^ = -5ew \C'\ e'^'^^'^^ , (21) 

where both the color-allowed and color-suppressed EWP amplitudes have approximately the 
same proportionality constant 

(5w ^ 0.65 ±0.15 . (22) 

These relations determine both the magnitudes and phases of the EW penguins. Their weak 
phases are equal to the weak phase of P', i.e. to — tt. They appear as the minus signs in 
Eqs. ipUl) and (j^H)- We do not use Eq. as a constraint in our fit, but simply use 6ew 
as a fit parameter and check whether it comes out within the expected bounds. 

Eqs. fl20|) and ()21|) were incorporated into Fit II but Fit I only employs redefined T' and 
C that cannot be directly related to the EW penguins. Instead, we write 

P'ew + P'ew = '^ew {T + C) e"*^ = -5ew {T + &) e"*^ (23) 
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and then neglect P'ew which is expected to be the smaller of the two to obtain 

~ -6ew [f' + C') e-'"' = ~6ew (|T'|e*^^ + \C'\e'^'^^'^^) . (24) 



This relation for P'ew used in Fit I while P'ew set to zero. 

The fitting parameters of both fits are shown in the columns for Fit I and Fit II in Tab 
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IVl An unusually large \C /T\ ^ 1.4 ratio predicted by Fit I is an indication of large \Ptu 
destructive interference between T and Ptu contributions to the redefined tree amplitude T, 
and constructive interference between C and Ptu contributions to C. Indeed, Fit II which 
separates Ptu and tree- level amplitudes predicts \Ptu\ = 14.9 and a much more reasonable 
\C/T\ = 0.46 ratio. 

Fits I and II represent a completely satisfactory description oi B ^ tttt and B — >■ nK 
decay modes. The branching ratio for B^ ir^K^ is predicted to be about 1.7 a below the 
observed value, while that for B^ 7i^K~^ is predicted to be about l.lo" below experiment. 
These deviations could be hints of new physics 0| , or simply due to underestimates 
of neutral-pion detection efficiencies 0l- The predictions are shown in the columns for 
Fits I and II in Tables IVII and I VI II Uncertainties for all predictions have been estimated 
by scanning the parameter space and studying the parameter sets that led to values no 
more than 1 unit above the minimum. The spread in predictions corresponding to those 
parameter sets has determined the uncertainties in predictions. The same method was used 
in an earlier analysis of i? ^ VP decays P]- 

The confidence level of Fit II is slightly lower than in Fit I because two new parameters 
{\Ptu\ and its strong phase 6p^^) have been added without a corresponding improvement in 
the value. The dependence of the weak phase 7 in Fits I and II is shown as the 

dotted and dash-dotted curves, respectively, in Fig. [TJ 
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TABLE V: Comparison of parameters extracted in fits to branching ratios and CP asymmetries 
under various assumptions. Values of the topological amplitudes arc quoted in units of cV. The 
Fit I column shows values for T and C and their strong phases in place of T and C amplitudes 
and phases. Probabilities are those for to exceed the value shown for the indicated number of 
degrees of freedom. 



Quantity 


Fit to 


7r7r,7ri^ 


Global fit 




Fit I 


Fit II 


TTTj- TTT 

J:* it ill 


J:* it iV 


7 


(61+27)° 


(65+35)° 


(66I16) 


/ cr y< +18\o 

(54I24) 


ITI 
1 1 


16.1+?° 

— i.y 


30.4+^V 


27.5 lb 3.2 


J +7.9 

27.4I4;g 




(34+25)° 


(17+?^)° 

\ — Iz / 


(25 ± 9)° 


(34l}D 


\c\ 


22.9+ti 

''—3.4 


13.9+^-0 

— 8.5 


19.2^|i 


24.3+5;^ 


Sr 




(-94+1^)° 


(-94I11) 


(-103^21) 


\P'\ 


48.2+°-9 


47.7+[]-^ 


4:7.7 ± U.9 


/IT Q+0.9 

47.8_i,i 


\Ptu\ 


(input) 


14.9+^^7° 


11.2 ± 3.4 


i2.3t^:^ 




(input) 


(3+9?)° 


(21 ± 16)° 


(37li^)° 


\s'\ 


(input) 


(input) 


09 1+3.0 

"^^•-L-3.3 


32.4I3.2 


Ss 


(input) 


(input) 


(-69_8 ) 


(' Tn+10\o 
(-70_8 ) 


\Stu\ 


(input) 


(input) 


u (^input ) 


c T+5.5 
0./_4 1 


^Stu 


(input) 


(input) 


(input) 


(-61^1)° 




^•^^-0.33 


n 42+0-50 
U-^^-0.29 


n 47+0-32 
-0.30 




Fit properties: 










xVd.o.f. 


7.34/8 


6.97/6 


18.06/15 


15.95/13 


CL (%) 


50 


32 


26 


25 


Derived quantities: 










1 P' 1 




o f:+3.6 


3.6li:| 




Ip/C 1 


(input) 


1 0+3.1 
J^-"J-1.0 


2.ll}i 


3.4ii:i 


\C/T\ 


1 IO+0.40 
-'- ■-■•-•-(), :-!! 


10+"-+^ 


0.70 ±0.16 


0.89 ±0.21 
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25 50 75 100 125 150 175 

■y (degrees) 



FIG. 1: (x^)min, obtained by minimizing over all remaining fit parameters, as a function of the 
weak phase 7. Dotted curve: Fit I; dash-dotted curve: Fit II; dashed curve: Fit III; sohd curve: 
Fit IV. Vertical dashed lines show the boundaries of the favored 95% confidence level range of 7 

Q 

(39° — 80°) from fits to CKM parameters |3| based on other measurements. 
VI. INCLUSION OF MODES WITH t] AND t]' 

To enlarge the fit and discussion to decays involving rj or 77' in the final state, we include 
an additional singlet amplitude. It is represented by 

S' = -\S'\e'^^ , (25) 

which gives two more fitting parameters. The relation between S and 5" is the same as the 
one between P and P': 

The importance of the S' amplitude has been discussed in Refs. jll|, |l5|, |21| mainly to 



account for the large branching ratios of the rj'K modes. Moreover, we include the parameter 
Ptu and its associated strong phase 6p^^. The penguin contribution Ptu is apparently required 
by our fits to decay modes involving i] and i]'. For instance, in n'^r]^'^ the Ptu 

contribution is of the same order as the other terms and cannot be neglected. The results 
under these assumptions are given in the column for Fit III in Table The dependence 
on 7 is shown as the dashed curve in Fig. ^ 
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TABLE VI: Comparison of predicted and experimental branching ratios in units of 10^^ and CP 
asymmetries for AS = B —>■ PP decays. The predictions of Fits I and II for r] and r]' modes 
are not reliable and are given for comparison purposes only. CP asymmetries, when predicted, 
are displayed on second line for a decay mode, while asymmetries in curly brackets (when shown) 
correspond to S (second line) and A (third line). 





JT It to 




Global fit 


J-JA-Lld llilClitCli 




Fit I 


Fit II 


Fit III 


Fit IV 


average 




c 1 9+0.38 


c 1 1 +0.22 
'^•-'--'^-0.14 


r. 1 1 +0.33 
"^•^^-0.37 


r 1 q+0.23 
^■-•^^-0.22 


5.2 ±0.8 




-0.00 ± 0.00 


-0.00 ± 0.00 


-0.00 ±0.00 


-0.01 ±0.00 


-0.07 ±0.14 


K+RO 


1.14 ±0.04 






1 Q1+0.99 
-^••5^-0.36 


< 2.5 










4.58t°;i? 


4.9 ±0.9 




-0 07+°-°^ 


-0.40l°|? 


-0.39l°:l? 


-0.401H1 


-0.44 ±0.18 




q qcr+0.60 


0.84tg:?^ 


^•^^-0.31 


2 qcc+0.89 


2.4 ± 1.1 (< 4.5) 




-0 07+°-°^ 


-0.4ll°|3 


-0.10 ±0.10 


^•^"^-0.34 





5° 



TT ' TT 



4.5^ 



hO.23 

-0.28 



4.55 



+0.07 

0.06 



4.58 



+0.10 
-0.12 



4.58" 



-0.08 

-0.11 



_0 7Q+0-25 
'-'■'^-0.16 



0.34 



+0.02 
-0.07 



1 qcc+0.17 



_n 74+0.26 

'-'■'^-0.21 

0.33 ± 0.02 
1 04+0.10 

-■-•^^-0.18 



0.44t?; 



35 
02 



0.57^?; 



25 
30 



7r% 



7777 



77' 77' 



^2+°-°'^ 

^•^^-0.20 

1.06 ± 0.04 
0.69 ± 0.02 

"•''^-0.03 



0.53tH^ 

i.78i5:i 

1 iq+3.40 
-'-•-'-^-0.83 

^•"^'-0.39 



-0.7410:22 
0.31 ± 0.06 

i.97ig:i 

0.54t0:22 
0.56lH^ 



-0 sq+o-^'^ 

^■°^-0.06 



0.30 



+0.02 

-0.04 



1 Q7+0.14 
^^•^'-0.19 

' n J2+0-53 
'-'•^"'-0.83 

^•^^-0.24 



4.6 ± 0.4 
-0.70 ±0.30 
0.42 ± 0.19 

1.9 ±0.5 



1.29 



+0.42 
-0.32 



1.21 



+0.92 
-0.33 



1 1^+0.30 
J^-J^^-0.33 

1.34 ±0.18 



1.67 
1.59 

0..38 



+0.86 
-0.44 

-0.80 
-0.42 

-0.18 
-0.10 



0.68 



+2.66 
-0.47 



1.54 



+0.40 
-0.29 



0.6612:1^ 



9 CI +0.51 
^•^-■--0.36 



0.16 



+0.61 

-0.12 



0.97 



+0.16 

-0.11 



0.95l°:?i 

9 1 f;+0.87 

0.68 ± 0.32 



< 1.5 

< 2.5 

< 3.7 

< 2.8 

< 4.6 

< 10 
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TABLE VII: Comparison of predicted and experimental branching ratios in units of 10^^ and CP 
asymmetries for \AS\ = 1 B —>■ PP decays. The predictions of Fits I and II for r] and rj' modes 
are not reUable and are given for comparison purposes only. CP asymmetries, when predicted, 
are displayed on second line for a decay mode, while asymmetries in curly brackets (when shown) 
correspond to S (second line) and A (third line). 



Mode 


Fit to TTTT, irK 


Global fit 


Experimental 




Fit I 


Fit II 


Fit III 


Fit IV 


average 






22.64lO:i 


22.05l°:i 


22.30t°:?^ 


21.8 ± 1.4 







0.00 ± 0.04 


1 n no 
'-'•^•J-0.03 


1 n no 
^•^^-0.03 


0.02 ± 0.06 




ii.+U_o.7o 


11 40+0-27 


11.40t°;?? 


li.OO„Q gg 


12.5 ± 1.0 




U-'J^-0.04 




0.07 ± 0.02 


0.091H3 


0.00 ±0.12 


r}K+ 


^•J^"-0.09 


o.2iini 


44+0.60 


3.63 ± 0.59 


3.7 ±0.7 







-o.iol°:°« 


-0.411H1 


-0.341°:^^ 


-0.52 ± 0.24 


rj'K+ 


29.38l?|? 


30.72^^:°^ 


74.56lJ;i 


75.2llJ;^^ 


77.6 ± 4.6 




0.01 ±0.01 


m+O-04 


0.02 ± 0.01 


0.01 ±0.02 


0.02 ± 0.04 






18.60l°i? 


18.89l°:ll 


18.78l°:^« 


18.2 ±0.8 




-0.101°:°? 


-0.10 ±0.01 


-0.10 ±0.02 


-o.ioin? 


-0.09 ± 0.03 




q 00+0.67 


9.29l°:io^ 


^•^^-0.65 


q 00+0.64 
^"''^-0.63 


11.7 ± 1.4 


< 


0.83 ± 0.01 

< 

i-o.iilHt 1 


' 0.831°:°^ J 
< 

-0 11+0-06 


0.83 ±0.01 

< 

[-0.121°:°^ 1 


0.83 ± 0.01 

< 

-0 11+0-05 


0.48 ± 0.42 
-0.40 ±0.29 




0.071H? 




2.66l°;« 


2.49lg;^f 2.5 ± 1.0 (< 5.2) 



-0.59l°:S 

o.6otO:i 

27.931?:^^ 

70+0-01 
^- 'U-0.00 



o.34i°:i 

, 0.85l°i| 
29.88tJ:|^ 



0.53lH^ 
0.021H1 



69.29 



+1.45 

-1.84 



0.04 



-0.00 
-0.02 



0.81: 



h0.09 
-0.05 



0.04 ± 0.06 



0.74 ± 0.01 

0.07 ± 0.02 



0.56ini 
, 0.03lHi 
69.271^:^^ 

n 7c;+0.00 



0.06 



+0.01 

-0.02 



65.2 ±6.2 
0.27 ±0.21 
-0.04 ± 0.13 
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Finally, since there is no reason to exclude such a term, we include a contribution from 
a singlet-penguin amplitude Stu associated with intermediate t and u quarks, consisting of 
a parameter \Stu\ and its associated strong phase 6stJ 

Stu = |5i,|e^('5s-+^) , (27) 

This exercise is denoted by Fit IV. The sole improvement with respect to Fit III is a better fit 
to the 71^1]' branching ratio, as shown in Table IVTl The tree amplitude |T| extracted 

from both Fit III and Fit IV is in agreement with the estimate obtained from a recent 
application of factorization [54] to the spectrum in i? ^ ttIu [s^, which yields 24.4 ±3.8 eV. 

Both Fit III and Fit IV represent a good description oi B ^ PP decay modes, including 
those with rj or r]' in the final state. The only problematic data points are the branching 
ratio for B^ n^K^ which is predicted to be about 1.7cr below the observed value and 
the mixing-induced asymmetry S{ri'K^) with the prediction (~ sin 2/3) at about 2.2a above 
the experimental value. The predictions for all other observed t] and t]' modes reproduce 
experimental values within their uncertainties. The predictions for as-yet-unseen modes are 
consistent with the current experimental upper limits on their branching ratios. The predic- 
tions are shown in the columns for Fits III and IV in Tables IVII and IVIII The dependence 
of on 7 in Fit IV is shown as the solid curve in Fig. ^ 

VII. STABLE AND LESS-STABLE ASPECTS OF FIT 
A. Robust aspects 

The value of the weak phase 7 obtained in B ^ PP data is consistent with other 
determinations. All versions of the fits have a local minimum in the range 48° < 7 < 73° 
(68% c.l.) allowed by global fits to phases of the CKM matrix |3| and near the range (63±6)° 
obtained in a fit to 5 —> VP data . The variation of 7 from fit to fit is at most about 12 
degrees, providing some idea of the systematic error associated with this approach. 

All fits are comfortable with a relatively large negative value of S'tttt which is the average 



of the Babar 



and Belle 



4l| values. Large negative 5',^^ is associated with larger a and 
smaller 7 (see, e.g., the plots in Refs. jsfi^). 



17 



The magnitude \P'\ of the strangeness-changing penguin amphtude changes very httle 
from fit to fit. ft is specified by the decay n^K^, which is expected to receive no 

other significant contributions. The presence of any direct CP violation in this decay would 
call that assumption into question, but no such asymmetry has yet been detected. 

All fits obtain a much larger value of \C/T\ than the range of 0.08 to 0.37 assumed in Ref. 



Moreover, all fits (including those to vrvr and ttK modes alone) entail a large strong 
relative phase 6c between the C and T amplitudes. The presence of a large color-suppressed 
amplitude is somewhat of a surprise from the standpoint of a priori calculations such as those 
in the QCD factorization approach and probably indicates a greater-than-anticipated 
role for final-state rescattering, which can generate such an effective amplitude (see also 
js?!). Such rescattering may be the reason why the decay —>■ tt^tt^ is more prominent 
than had been expected. All our fits now entail a branching ratio for this mode of about 
2 X 10~^. Although the favored values of some topological amplitudes (e.g., C, Ptu) show 
noticeable variations from fit to fit, they change together in a correlated way so that the 
predictions for almost all of the modes that involve them remain very stable. 

The strong phase 5t of the tree amplitude T with respect to the penguin amplitude P is 
found to be non-zero and of the order of 20° to 30°. It is most likely driven by the need to 
simultaneously describe a large direct CP asymmetry (the parameter At^.^) in B^ tt^tt' 
and a small but significant direct asymmetry in B^ 7i~K~^. These quantities are well-fitted 
and their predicted values do not differ much among the four fits. 

While the electroweak penguin parameter 6ew was initially constrained to lie within 
the range ()22|1 . we found that leaving it as a free parameter led to results consistent with 
that range except in the cases of Fit II and Fit III. Thus, our fits do not favor a large 
phenomenological EWP amplitude. This should be contrasted with Refs. 0, where a 
different assignment of weak and strong phases is given in expectation of new physics con- 
tributions. Our fits also do not favor much deviation of the predicted S^j^oko time-dependent 
asymmetry parameter from its predicted standard-model value of sin(2/3) ~ 0.74 59l |. 

Once one admits enough parameters into the fits to correctly describe modes involving 
rj and rj', the negative direct CP asymmetry in -K^rj observed by BaBar is 

correctly reproduced. The possibility that this asymmetry could be large was first noted in 
Ref. j60| and pursued in Refs. 15|. We predict a similarly large negative CP asymmetry 
in B^ — > T]K~^, as observed 30|. These asymmetries can be large because no single weak 
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amplitude dominates the decays. As sensitivities of asymmetric e^e~ collider experiments 
improve through the accumulation of larger data samples, we expect more such decay modes 
to emerge. 

The mixing-induced and direct asymmetries S{ri'K^) and A{r]'K^) are predicted to be 
close to sin(2/5) and 0, respectively. These two values would be expected if the — >• 
rj'K^ decay amplitude had consisted of just QCD penguin P' and singlet penguin S'. The 
interference of these terms with the much smaller C, Pj.^, and S'^^ amplitudes leads to small 
deviations from the expected values. These deviations are to a large extent determined by 
the ratio \Aq/A'p\ of the terms with the weak factor V*f^Vus (C, P/„, and S'{„) and the terms 
with the weak factor V*ijVcs {P' and S'). ^An/A'A is typically predicted by QCD factorization 



and PQCD to be smaller than 0.02 



12,|6lll63]- Our best conservative estimate of \A(j/Ap\ 



is based on Fit IV. We find that the SU(3) fit prefers somewhat larger values: \A'(j/ A'p\ = 
OM2toml Fit HI (somewhat more stable than Fit IV) predicts \A'c/A'p\ = OMOtoml 
More conservative bounds on \A'^/A'p\ and on the asymmetries Siv'K^) and A{ri'K^) were 
obtained recently in a model-independent way using flavor SU(3) (63l |. 

We have explored the effects of changing the rj-r]' octet-singlet mixing angle from its 
nominal value 6 ~ 19.5° defined in Sec. II. The angle 6 assumed a value of 22.0° in Fit III 
with a free mixing angle while of the fit improved by just 1.12. With one additional 
parameter in the fit, this did not result in a better fit quality. Fit IV with a free mixing 
angle preferred 6 = 20.4°, with the fit quality dropping by 5%. Thus, leaving the rj-r]' 
mixing angle as a free parameter, we found variations of only a few degrees and negligible 
improvements in fits. 



B. Aspects sensitive to assumptions 

The possibility of a large Ptu term in Fit II leads to a wide range of predicted branching 
ratios for — ^ K^K^ and B^ K^K^. This range is considerably reduced in other fits. 

The magnitude and phase of the singlet penguin amplitude S' are probably not well- 
determined. The two quantities are correlated, as first pointed out in Ref. ^] and noted 
further in Ref. (ll|. For example, a much smaller magnitude of S' is required to fit the 
charged and neutral B — > rj'K decay modes if 5" and P' (the gluonic penguin amplitude) in- 
terfere constructively with one another. The QCD factorization approach |lj] finds negligible 
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S' contribution to these decays, explaining their enhancement by means of nonet-symmetry 
breaking effects as proposed, for example, in Ref. 6^, and making use of the constructive 

interference of non-strange and strange components of the rj' in the gluonic penguin ampli- 

|i i 

tude |65']. One should also point out that many other explanations have been proposed for 
the enhancement of -B ^ rj'K modes Ig^. One also finds the magnitude of S' to be 
sensitive to small changes in the octet-singlet mixing in rj and t]'. 

Predictions for the branching ratio and CP asymmetry in rjK^ depend crucially 

on the introduction of the S' amplitude. Since this amplitude is uncertain in magnitude 
and phase, those predictions (although apparently satisfied) should be viewed with caution. 
The same warning applies to the mode — > rj-K^ . 

As already noted, the predicted branching ratio for B^ 71^1]' is quite sensitive to 
assumptions, and was the sole quantity which could be compared to experiment that led to 



the introduction of the Stu term in Fit IV. In Ref. ll| we noted a tight correlation between 
predicted branching ratios and CP asymmetries for B^ — > ii^r] and B^ -n^rj . With the 
added possibility of nonzero Ptu and Stu contributions, this correlation no longer holds. 

The only other prediction whose values are significantly different in Fits III and IV is 
the mixing-induced asymmetry 5'(7r*^7r'^). One should trust the larger values of this quantity 
predicted by Fits I and II. These fits to vrvr and t^K data points are not affected by the 
uncertainties associated with r/ and Vj . Their predictions for the asymmetries in 5° 



0^0 



7r TT 

modes thus are expected to be more rehable. 

The introduction of the Stu term changes the favored value of 7 by a noticeable amount, 
though still within limits from CKM global fits j^. As noted, this provides one estimate of 
systematic errors associated with analyses of the present form. 

VIII. MODES TO BE SEEN 

Several decay modes are predicted to occur at levels just below present upper bounds, 
and can provide useful constraints on the residual uncertainties in our fits. For example, 
the decays B^ —>■ K^T^ and B^ K^lC^ are predicted to have branching ratios exceeding 
10^^ (somewhat larger than in Ref. with the exact value depending on the fit. The 

decay modes B^ — > tt^i] and B^ — > Tc^rj' also should be visible at this level. The modes B^ — > 
[riT], T]-)]', t]'-)]') will probably require more work. We also make predictions for the direct and 
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TABLE VIII: Comparison of observed and predicted direct CP asymmetries for some B — > tttt 
and B ttK decay modes. 



Decay Exptl. Present QCDF [12j PQCD 

Mode average work (a) Full range Favored (b) 



50 



B' 



0.42 ±0.19 0.30: 



-0.02 

-0.04 



-0.07; 



-0.14 
-0.13 



B+^-K-K+ -0.09 ±0.03 -0.10l[!:J]^ 



-0.40 ±0.29 -0.11 



4-0.05 
-0.02 



04+°-0^ 

'-'•'-'^-0.10 



-0.03 ±0.04 



0.10 
-0.04 
0.01 



0.16 to 0.30 
-(0.13 to 0.22) 



Ref. \YJj 



Input 
-n 14+0.09 



-0.05 



-fO.29 
-0.24 



(a) Fit IV; (b) Scenario S4 



TABLE IX: Comparison of observed and predicted direct CP asymmetries for some B decay modes 
involving r] and r/'. 



Decay 
Mode 


Exptl. 

average 


Present 
work (a) 


QCDF \12\ 

Full range 


Favored (b) 


B+ TT+n 
B+ r]K+ 


-0.44 ±0.18 
-0.52 ±0.24 


-0.40t°;°^ 
-0.34^°;J^ 


-0.15 ±0.20 
-0 1Q+0-29 

^■-'■^-0.30 


0.06 
0.10 



(a) Fit IV; (b) Scenario S4 

mixing-induced asymmetries in -k^-k^ and — > rjK^ , with yl(7r°7r°) exceeding 0.5. A 

prediction for the branching ratio of K^K^ cannot be made in our approach. The amplitude 
of this decay mode receives contributions from exchange and penguin annihilation diagrams 
that are neglected in this paper. It is very desirable that a more strict experimental upper 
limit be set for this mode to justify the assumption of negligibility of similar contributions 
to other neutral AS* = decay modes. 



IX. COMPARISON WITH OTHER APPROACHES 

The signs of our predicted direct CP asymmetries agree with those measured experimen- 
tally for the five processes in which non-zero asymmetries are reported at greater than the 
2(T level. We summarize these and our predictions for them in Tables IVIIII and IIXI (For 
others, as shown in Tables VI and VII, negligible asymmetries are predicted, in accord with 
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observation.) 

The fact that we agree with all five signs and magnitudes is due in part to the fiexibility 
of our SU(3) fit, but still represents a non-trivial consistency in our description of strong 



phases. We were not able to achieve this consistency in Ref. The same correlation 

between predicted signs of direct asymmetries in tt+tt^ and 5° n^K^ occurs in all 

the methods compared in Table IVIIII A definite prediction of the absolute signs, in accord 
with experiment, is made in Ref. 

Fits to B ^ PP branching ratios in the various approaches which we compare with ours 
[l^ llfil . IitI lisj l are generally acceptable, especially when allowance is made for possible large 
penguin amplitudes and color-suppressed contributions. These fits now are converging on 
a preference for 7 in the range preferred by fits to other observables constraining CKM 
parameters. 



X. SUMMARY 



We have analyzed the decays of B mesons to a pair of charmless pseudoscalar mesons 
within a framework of fiavor SU(3). Acceptable fits to B ^ titi and B Ktt branching 
ratios and CP asymmetries were obtained with tree, color-suppressed (C), penguin (P), 
and electroweak penguin (Pew) amplitudes, but in order to describe processes involving 
7] and 7]' we needed to include a large fiavor-singlet penguin amplitude (5) and a penguin 
amplitude Ptu associated with intermediate t and u quarks. For the B^ -n^rj mode a 
term Sm associated with intermediate t and u quarks also was employed. 

We were able to achieve a good fit to the five most significant direct CP asymmetries, 
as noted in Tables IVlllI and IIXI We found values of the weak phase 7 roughly consistent 
with those obtained earlier in an analysis of -B — > VP decays [(7 = 63 ±6)°], and with other 
analyses j^l of CKM parameters, for which the 68% confidence level limit is 48° < 7 < 73°. 
A global fit without Stu gave 7 = (66li6)°, while adding Stu yielded 7 = The 
difference between these two serves as an estimate of systematic error. 
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